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SUNMARY 


This  report  describes  an  experimental  and  analytical  evaluation  of  the 
vibrational  characteristics  of  a  truncated  conical  shell  of  sandwich 
construction.  The  shell  was  composed  of  fiber  glass-epoxy  facings  and 
aluminum  honeycomb  core.  The  nominal  shell  geometry  had  a  5-degree  semi¬ 
vertex  angle,  a  diameter  of  58  inches  at  the  large  end,  and  a  length 
of  72  inches.  These  dimensions  were  selected  to  be  typical  of  an  Army 
aircraft  fuselage  structure.  The  shell  was  suspended  by  soft  cords,  to 
closely  approximate  free-free  boundary  conditions,  and  was  excited  laterally 
by  an  electrodynamic  shaker  placed  at  various  locations  along  the  shell. 

The  shell  surface  was  instrunented  with  600  metallic- foil  strain  gages. 
Accurate  separation  of  vibrational  modes  having  closely  spaced  frequencies 
was  accomplished  by  means  of  the  Kennedy-Pancu  method  of  data  reduction, 
using  measured  frequency,  strain  amplitude,  and  relative  phase  angle 
between  response  and  excitation.  The  two  lowest  unsymmetric-mode  fre¬ 
quencies,  for  various  values  of  circumferential  wave  number,  agreed  quite 
closely  with  those  predicted  by  the  Rayleigh-Ritz  inextensional  analysis 
presented  in  the  report.  This  analysis  includes  the  orthotropic  nature 
of  the  compos ite-material  facings.  Values  of  the  damping  logarithmic 
decrement  are  also  presented. 
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This  report  was  prepared  by  the  University  of  Oklahoma  Research  Institute 
(OURT)  under  Phase  I  of  U.S.  Army  Aviation  Materiel  Laboratories  (USAAVLABS) 
Contract  DAAJ02-67-C-0111  for  research  conducted  during  the  period  from  July 
1,  1967  to  July  31,  1968.  Work  under  this  contract  has  also  included  bi¬ 
axial  loading  of  fiber- reinforced-plastic  laminates  (Phase  II),  reported 
separately  in  Technical  Report  68-86,  titled  Behavior  of  Fiber-Reinforced- 
Plastic  Laminates  Under  Uniaxial,  Biaxial,  and  Shear  Loadings. 

The  research  effort  is  a  continuation  of  the  research  accomplished  under 
three  previous  contracts  and  reported  in  the  following  USAAVLABS  reports: 

1.  Technical  Report  67-65,  Fabrication  and  Full-Scale  Structural 
Evaluation  of  Sandwich  Shells  of  Revolution  Composed  of  Fiber 
Glass -Reinforced  Plastic  Facings  and  Honeycomb  Cores. 

2.  Technical  Report  65-66,  The  Effect  of  Resin  Content  and  Voids  on 
the  Strength  of  Fiberglass -Reinforced  Plastics  for  Airframe  Use. 

3.  Technical  Report  65-60,  Dpamic  Elastic,  Damping,  and  Fatigue 
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Reinforced  Sandwich  Structure  for  Airframe  Use. 

The  present  report  was  written  by  Dr.  Charles  W.  Bert,  project  director  and 
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1.  INTRODUCTION 


Under  the  sponsorship  of  the  U.S.  Army  Aviation  Materiel  Laboratories 
(USAAVLABS) ,  the  University  of  Oklahoma  Research  Institute  (OURI)  has 
conducted  a  major  research  program  for  the  past  five  years  in  the  area 
of  fiber-glass  reinforced  plastics  (FRP)  suitable  for  primary  structures 
of  U.S.  Army  aircraft.  Tasks  which  have  been  completed  include  the 
evaluation  of: 

(1)  Several  sandwich  fabrication  techniques  (Reference  1) . 

(2)  Effect  of  primary  fabrication  variables  on  conventional  strength 
properties  of  laminates  (References  1,  2,  3). 

(3)  Several  honeycomb- core  sandwich  configurations  as  small  panels 
sized  to  fail  in  face  rupture  and  various  buckling  modes  (References 
1,4,5). 

(4)  Fatigue  properties  of  sandwich  material  as  beams  with  FRP  facings  and 
honeycomb  core  (References  6  and  7). 

(5)  Dynamic  stiffness  and  damping  of  sandwich  material  in  the  form  of 
beams  (References  6  and  8). 

(6)  Fabrication  and  static  buckling  testing  of  full-size  sandwich  shell 
structures  in  the  foim  of  cylindrically  curved  panels,  complete 
cylinders,  and  truncated  cones  (References  9,  10,  11,  and  12). 

(7)  Mechanical  behavior  of  flat  laminates  under  uniaxial  tension,  biaxial 
tension  (1:1  and  1:2  principal- stress  ratios),  and  pure  shear 
(Reference  13). 

The  objectives  of  the  present  research  were  (1)  to  measure  resonant  fre¬ 
quencies  and  associated  modal  strain  distributions  and  damping  factors  of 
a  truncated  conical  sandwich  shell  suspended  in  the  free-free  condition 
and  (2)  to  compare  these  quantities  with  theoretical  predictions. 

A  nunber  of  experimental  investigations  have  been  carried  out  to  evaluate 
theories  of  vibration  of  sandwich- type  beams  and  plates  (Reference  14). 
However,  the  authors  do  not  know  of  any  investigation  on  sandwich-type 
shell  structures,  which  are  widely  used  in  aircraft  structures.  The  shell 
configuration  used  in  the  present  investigation  was  that  of  a  truncated 
cone.  The  scale  selected  was  sufficiently  large  to  be  typical  of  an  Army 
aircraft  fuselage  and  to  achieve  lowest  resonant  frequencies  sufficiently 
low  to  enable  accurate  measurement  of  a  large  nunber  of  natural  modes. 

Also,  it  is  difficult  to  scale  down  a  sandwich- type  structure;  i.e.,  if  the 
core  thickness  is  too  thin,  the  sandwich  effect  is  negligible. 

In  the  interest  of  economy,  the  shell  was  fabricated  with  the  same  equip¬ 
ment  and  techniques  previously  developed  by  OURI  for  the  shell  buckling 
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experiments  (Reference  9).  The  inside  and  outside  facings  consisted  of 
two  plies  of  828- Z  epoxy  resin  resin  reinforced  by  Vo7 an-A- finished  E- 
glass  in  the  form  of  181-style  fiber  glass  fibric.  The  core  material  was 
5052  aluminum- alley ,  hexagonal-cell  nonperforated  honeycomb. 

To  be  auite  certain  that  the  boundary  conditions  achieved  in  the  experiments 
matched  those  used  in  analysis,  the  boundary  conditions  selected  for  the 
experiments  were  free  edges.  Also,  these  conditions  were  the  least  expen¬ 
sive  to  achieve  experimentally  and  were  used  in  previous  experiments  on 
homogeneous  shells  at  the  Langley  Research  Center  (References  15,  16,  and 
17)  and  at  Southwest  Research  Institute  (Reference  18). 

Although  the  free-free  condition  is  most  easily  achieved  experimentally,  it 
is  perhaps  the  most  difficult  condition  frem  an  analytical  viewpoint, 
particularly  in  the  case  of  a  sandwich  shell,  because  of  the  difficulty  in 
finding  functions  tliat  satisfy  the  boundary  conditions.  This  is  attested 
to  by  the  very  limited  nunber  of  analyses  of  homogeneous  conical  shells 
(or  even  homogeneous  cylindrical  shells)  with  these  boundary  conditions. 

It  can  be  applied  to  truncated  conical  shells  of  either  simple  single- layer 
or  sandwich  construction  with  either  isotropic  or  orthotropic  (composite- 
material)  facings.  Inextensional  motion,  in  which  the  extensional  strains 
are  neglected,  is  assumed.  Solution  is  carried  out  by  the  Rayleigh-Ritz 
energy  method  for  the  twe  inextensional  modes,  which  are  the  lowest  unsym- 
metric  modes.  The  circumferential  wave  number,  n,  may  be  any  integer 
greater  than  one.  The  validity  of  the  simplifying  assumptions  made  in  the 
analysis  are  verified  by  the  good  agreement  between  the  experimental  and 
analytical  results. 
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2. 


SPECIMENS  AND  MATERIAL 


Two  large-scale  sandwich  conical-shell  specimens  were  used  in  the  experi¬ 
mental  program.  The  first  specimen,  referred  to  as  the  preliminary  shell, 
was  used  for  instrumentation  and  test  procedure  checkout  only.  This 
specimen  was  made  by  cutting  off  the  locally  damaged  region  from  the  last 
shell  tested  in  the  sandwich  shell  buckling  program  (Reference  9) ,  thus 
making  it  approximately  6.5  inches  shorter  than  the  final  shell. 

The  second  specimen,  referred  to  hereafter  as  simply  the  shell,  was  made 
on  the  same  mold  as  the  other  conical  shells  (Reference  9).  However,  in 
order  to  be  able  to  handle  the  prepreg  more  easily  during  fabrication,  the 
B-stage  time  was  extended  to  14  hours  for  the  final  specimen.  This  was 
double  the  B-stage  time  required  to  obtain  optimal  strength  properties,  but 
high-strength  properties  were  not  required  in  the  present  program.  The 
resulting  shell  had  considerably  fewer  wrinkles  than  those  used  in  the 
buckling  program. 

Table  I  gives  data  on  the  sandwich  consitutents,  and  Table  II  lists  the 
dimensional  and  weight  characteristics  of  the  shell.  The  material  proper¬ 
ties  of  the  facings  and  core  are  listed  in  Table  III.  The  facing  tensile 
properties  and  shear  modulus  were  determined  from  tests  on  small  tensile 
and  torsion-iube  specimens,  which  were  subject  to  the  same  curing  cycle 
as  the  facing  laminates. 


TABLE  I. 

SANDWICH  CONSTITUENTS 

-  1 

Constituent 

Description 

Thickness  Orientation 

Facing 

82  8- Z  Epoxy, 

181-E  Volan  A 

Fiber  glass 

Two-ply 
(0.02  inch 
total) 

Warp  parallel  to 
shell  axis 

Core 

5052  Aluminum, 

1-mil  Nonperf orated 
Foil,  1/4- inch 
Hexagonal  Cell 

0.3  inch 

Ribbon  direction 
parallel  to  shell 
axis 

Adhesive 

AF-110B  Film- 
Supported  Epoxy 

0,015  inch 

- 

3 


TABLE  II.  FINAL  SANDWICH  SHELL  DIMENSIONS  AND  WEIGHT 


3.  EXPERIMENTAL  EQUIPMENT 

The  shells  were  suspended  from  a  large  steel  frame  by  six  soft  springs 
(Figure  1)  so  that  the  suspension  resonant  frequency  was  below  1  cps.  The 
frame  was  previously  used  for  casting  molds  (Reference  9).  An  electro¬ 
dynamic  exciter  (MB  Model  Cll,  50  pounds  maximum  force  capacity)  was  attached 
to  the  shell  by  a  force  link  that  was  designed  to  reduce  the  coupling  be¬ 
tween  the  specimen  and  the  steel  frame.  The  exciter  was  attached  at  differ¬ 
ent  locations  around  the  shell  and  along  the  length  to  determine  the  exciter 
effects  on  the  data. 

The  instrumentation  transducers  used  to  measure  the  strain  distribution  at 
the  resonant  frequencies  were  600  metallic- foil  strain  gages  (Budd  Model 
C6-141B)  located  on  the  outer  facing  of  the  shell  in  an  array  as  shown  on 
Figure  2.  Two  gages  were  located  at  each  grid  point,  one  to  measure 
circumferential  strain  and  one  to  measure  meridional  strain. 

The  electronic  circuitry  that  was  used  to  monitor  and  measure  the  dynamic 
signals  from  the  gages  is  shown  schematically  in  Figure  3,  and  a  photograph 
is  shown  in  Figure  1.  The  strain  gage  was  used  as  a  variable  resistor, 
whose  resistance  change  was  a  function  of  the  elongation  applied  to  the 
gage.  A  voltage  of  7  volts  D.C.  was  applied  to  the  gage  from  a  300-volt 
D.C.  power  supply  through  a  variable  lOk-ohm  dropping  resistor.  Analysis 
of  the  circuitry  shows  that  net  only  was  this  a  voltage-dropping  resistor, 
but  it  was  also  a  high  impedance  source  to  the  A.C.  signal  caning  from  the 
strain  gage.  The  signal  was  applied  to  the  input  of  two  A.C.  amplifiers 
(Eico,  Model  250)  in  series  which  had  a  high  impedance  to  the  D.C.  portion 
of  the  signal,  allowing  only  the  A.C.  component  to  be  amplified.  Fran  the 
output  of  the  amplifier,  the  signal  was  supplied  to  various  instruments. 

The  voltmeter  (Hewlett-Packard  Model  400H)  and  oscilloscope  (Tektronix 
Model  503)  were  used  to  monitor  the  strain-gage  amplitude  and  wave  shape. 

Each  strain-gage  system  was  calibrated  separately.  A  known  sinusoidal 
signal  was  applied  across  each  strain  gage  as  indicated  in  Figure  3.  With 
a  known  input  amplitude,  the  signal  was  recorded  on  the  oscillograph.  The 
300-volt  D.C.  power  supply  was  removed  when  the  systems  were  being  cali¬ 
brated.  With  the  system  calibrated  in  this  manner,  errors  in  the  absolute 
strain  values  were  minimized. 

The  frequency  response  of  the  amplifier  recording  system  has  a  flat  response 
from  10  to  25,000  cps;  the  range  of  this  experiment  was  from  5  to  500  cps. 

The  responses  of  the  strain  gages  for  dynamic  measurements  were  not  measured 
in  this  test,  but  according  to  Reference  19,  the  responses  are  from  0  to 
20,000  cps. 

The  equation  for  the  strain  measurements  taken  from  the  strain-gage  system 
is  (Reference  20) : 


strain  =  c  =  (/\R'/R')/(G.F.) 
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Figure  1.  Experimental  Vibration  Test  Setup  (A,  Conical  Sandwich  Shell; 
B,  Soft  Suspension  Springs;  C,  Vibration  Exciter;  D. 
Instrumentation  Cart) . 
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Figiirc  2.  Layout  of  Shell  Showing  Strain-Gage  Locations. 
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where  G.F.  =  gage  factor 

AR’  =  change  in  gage  resistance 
R'  =  gage  resistance. 

Analysis  of  the  electronic  circuitry  indicates  that 

(AR'/R')  =  (AV/V) 

where  AV  =  change  in  voltage  across  the  gage 

V  =  D.C.  voltage  across  the  gage. 

Therefore , 

ed  =  (AV/V) (G.F.) 
where  =  dynamic  strain. 

The  voltage  change  was  taken  from  the  strain-gage  output  as  noted  on  the 
voltmeter  after  applying  the  appropriate  calibration  factor.  Substitution 
into  the  above  equation  gave  the  djmamic  strain  values. 

An  accelerometer  (Endevco,  Model  FA-72)  was  attached  to  the  armature  of 
the  vibration  exciter  to  provide  a  measure  of  the  input  force  into  the  shell. 
The  output  of  the  accelerometer  (Endevco,  Model  2016)  was  amplified  and 
then  was  fed  into  a  readout  system,  as  indicated  on  Figure  3.  The  voltmeter 
was  an  indication  of  the  output  level,  and  the  oscilloscope  was  used  to 
monitor  the  signal  to  assure  that  there  was  no  distortion. 

The  accelerometer  system  was  calibrated  by  introducing  a  known  signal  in 
series  with  the  accelerometer  and  by  recording  the  output  level.  This  was 
done  in  accordance  with  Reference  21. 

The  outputs  of  the  strain- gage  system  and  the  accelerometer  were  applied  to 
each  side  of  the  digital  phase  meter  (AD-YU  Digital  Phase  Meter,  Model  524A 
and  EAI  Digital  Voltmeter,  Model  5002A) ,  which  indicated  the  relative  phase 
between  the  input  and  strain  distribution. 

The  system  indicated  on  Figure  3  is  for  the  output  amplitude  and  phase  of 
a  single  strain  gage.  Only  one  gage  could  be  monitored  at  a  time,  since 
only  one  data  acquisition  system  was  available.  A  telephone  switchboard 
was  used  to  change  gage  locations.  The  switchboard  was  found  to  give  the 
best  signal- to-noise  ratio  at  the  lowest  cost;  other  fas  ter- switching 
arrangements  were  tried,  but  each  increased  the  noise  level.  The  strain- 
gage  signals  were  on  the  order  of  100  microvolts,  which  made  the  noise 
level  of  prime  importance. 

Other  instruments  of  minor  importance  were  included  in  the  data  acquisition 
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system.  An  oscillograph  (Honeywell  Visicorder,  Model  906A)  was  used  to 
record  the  log-decrement  damping- factor  signal.  In  this  system,  the  signal 
was  paralleled  into  the  oscillograph  from  the  voltmeter. 

A  stroboscope  (General  Radio  Model  1531-AB)  was  used  sometimes  to  help 
define  the  modal  shape  of  the  shell  under  excitation. 

A  second  oscilloscope  was  used  to  display  a  Lissajous  pattern  as  an  indi¬ 
cation  of  the  phase  between  the  input  and  strain.  This  was  used  as  a 
check  on  the  phase  angle  quadrant  obtained  from  the  phase  meter. 

A  band-pass  filter  (SKL  Model  302)  was  used  sometimes  to  help  eliminate 
the  noise  in  the  strain-gage  signals.  The  filter  was  used  sparingly, 
since  it  introduced  a  phase  shift  that  was  dependent  on  frequency.  When 
the  phase  data  were  of  prime  importance,  the  filter  was  not  used. 
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4 .  EXPERIMENTAL  PROCEDURE 


Ihe  procedure  used  to  acquire  the  resonant- frequency  and  modal-strain- 
distribution  data  at  each  resonance  was  (1)  to  locate  the  resonant  fre¬ 
quency  using  a  modified  Kennedy-Pancu  technique  (see  the  Appendix  for 
development  of  this  technique)  and  (2)  to  define  the  modal  shape  at  this 
frequency. 

The  resonant  frequency  was  located  by  monitoring  a  selected  series  of  gages 
while  the  system  was  excited  throughout  the  frequency  range  of  interest. 

In  this  first  survey,  the  peak  amplitude  method  was  used  to  locate  the 
approximate  range  of  the  resonant  frequency.  The  phase  between  the  input 
and  gage  was  not  needed  in  this  survey.  The  gages  were  selected  by  exciting 
the  shell  and  by  noting  the  antinodal  points.  Inherently,  certain  portions 
of  the  shell  were  excited  at  all  the  resonant  frequencies.  The  survey  was 
made  by  monitoring  the  input  acceleration,  strain  outputs,  and  frequencies. 
Frequency  intervals  were  taken  at  5  cps  across  the  entire  band,  with  1  cps 
near  resonance.  Closer  intervals  were  taken  to  define  the  resonant  point 
completely.  Once  the  resonant  frequencies  were  located,  data  were  taken 
at  these  frequencies  to  make  the  modified  Kennedy-Pancu  plot. 

The  Kennedy-Pancu  data  were  taken  by  varying  the  excitation  frequency  and 
by  monitoring  the  phase  angle  between  the  input  acceleration  and  the  strain- 
gage  signal.  The  gage  selected  for  these  data  was  the  one  that  gave  the 
largest  output  in  the  preliminary  peak- amplitude  survey.  The  excitation 
frequency  was  varied  from  a  point  below  resonance  that  indicated  near- zero 
phase  to  a  point  above  resonance  that  indicated  that  the  phase  was  returning 
to  zero.  Since  the  phase  meter  would  not  function  below  a  threshold  input 
voltage  (0,35  volt),  data  were  taken  only  in  the  general  vicinity  of  reso¬ 
nances.  These  data  were  sufficient  to  define  the  uncoupled  resonances  by 
means  of  the  characteristic  circular  arcs  on  the  modified  Kennedy-Pancu 
plots.  To  minimize  the  nonlinear  effects  of  the  shell,  the  strain  level 
exhibited  on  the  strain  gage  was  held  constant,  and  the  input  level  into 
the  shell  was  varied.  The  excitation  frequency  was  changed  until  approxi¬ 
mately  10-degree  variation  in  phase  was  noted.  The  input  force  was  varied 
until  the  gages  exhibited  a  predetermined  output;  then,  the  input  acceler¬ 
ation,  output-gage  signal,  phase  relation,  and  frequency  were  recorded.* 

This  same  procedure  was  repeated  until  the  output  signal  became  too  small. 
The  modified  Kennedy-Pancu  plot  was  drawn  fran  these  data,  and  the  uncoupled 
resonant  frequencies  were  determined.  The  excitation  frequency  of  the  shell 
was  set  at  the  uncoupled  resonant  frequency,  and  the  excitation  level  was 
not  changed  from  the  previous  data.  Each  strain-gage  output  and  the  corres¬ 
ponding  phase  relation  were  read.  Fran  these  data,  the  modal  strain  distri¬ 
bution  was  determined.  Once  the  data  were  acquired,  the  shell  nodal  points 
were  investigated  and  checked  by  locating  the  points  by  feel  and  also  by 


*At  times,  the  phase  meter  would  not  give  an  accurate  reading  in  relation 
to  the  angle’s  quadrant.  The  Lissajous  pattern  was  used  as  a  check  to 
make  sure  that  the  quadrant  was  correct. 
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sight  using  the  stroboscope. 


The  same  procedure  was  repeated  for  each  resonant  frequency  throughout 
the  frequency  range.  Resonant  frequencies  were  investigated  from  5  to 
400  cps.  The  lower  limit  was  below  the  first  resonant  frequency  and  was 
the  lowest  frequency  of  the  exciter.  The  upper  limit  was  reached  when 
the  strain  signal-to-noise  ratio  was  too  low  to  obtain  reliable  data. 


5.  ANALYSIS 


a.  Review  of  Previous  Analyses 

Although  sandwich- type  construction  is  coming  into  extensive  use  in 
aircraft  structures  in  shell  configurations,  only  a  very  limited 
range  of  vibrational  problems  involving  such  structures  have  been 
solved.  Apparently,  the  first  such  analysis  was  performed  by  Yu 
(Reference  22)  for  the  free  vibrations  of  a  simply  supported  complete 
cylinder  of  sandwich  construction.  This  was  followed  by  work  on  the 
following  aspects  for  the  same  configuration:  wave  propagation  by 
Chu  (Reference  23),  large-amplitude  vibration  by  Chu  (Reference  24), 
harmonic  forced  vibration  by  Bieniek  and  Freudenthal  (Reference  25) , 
structural  damping  by  Yu  (Reference  26)  and  Jones  and  Salerno  (Refer¬ 
ence  27),  and  effect  of  initial  stress  by  Greenspon  (Reference  28). 

The  cylindrically  curved  panel  configuration,  rather  than  a  complete 
cylinder,  was  analyzed  by  Mead  and  Pretlove  (Reference  29)  and  by 
Jacobson  and  Wenner  (Reference  30).  Shallow  spherical- shell  caps  were 
analyzed  by  Tasi  (Reference  31)  and  by  Koplik  and  Yu  (References  32, 

33,  and  34).  Suvemev  (Reference  35,  as  reported  in  Reference  36) 
treated  conical-frustun  shells. 

None  of  the  analyses  mentioned  above  are  applicable  to  sandwich  shells 
with  orthotropic  facings,  encountered  in  practice  in  the  form  of  fiber- 
reinforced  composites.  However,  some  of  them  did  take  into  considera¬ 
tion  the  much  more  simple  effect  of  an  orthotropic  core,  as  characterized 
by  the  commonly  used  hexagonal- cell  honeycomb  core.  The  first  sandwich- 
shell  vibrational  analyses  to  consider  orthotropic  facings  were  made 
independently  by  Azar  (Reference  37),  Vasitsyna  (Reference  38,  as  re¬ 
ported  in  Reference  36) ,  and  Baker  and  Herrmann  (Reference  39) .  Azar 
treated  axisymmetric  free  vibrations  of  arbitrary  open-ended  shells  of 
revolution,  as  exemplified  by  conical  and  paraboloidal  shell  frusta, 
with  simply  supported  edges.  Vasitsyna  analyzed  free  vibrations  of 
simply  supported  circular  cylinders,  while  Baker  and  Herrmann  considered 
the  same  problem  with  the  addition  of  a  general  state  of  initial  stress. 
Later  Bacon  and  Bert  (Reference  40)  extended  Azar's  work  to  include 
unsymmetric  modes. 

It  is  important  to  mention  that  none  of  the  analyses  discussed  above 
considered  free-edge  boundary  conditions.  Most  of  than  used  simply 
supported  edges,  and  a  few  used  clamped  edges.  As  stated  in  the  intro¬ 
duction,  this  is  probably  due  to  the  increased  mathematical  complexity 
of  the  free-edge  condition.  In  fact,  only  a  very  limited  number  of 
analyses  of  homogeneous,  isotropic,  conical  shells  with  free  edges  have 
been  carried  out,  and  these  only  quite  recently:  Hu  (Reference  41); 

Hu  et  al  (Reference  42);  Sewall,  mentioned  in  a  report  by  Mixson 
(Reference  17) ;  and  Krause  (Reference  43) . 

Hu's  analysis  was  formulated  to  include  both  membrane  and  bending 
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effects,  with  solution  by  the  Galerkin  method.  However,  the  only 
numerical  results  which  he  published  for  free-free  boundary  conditions 
include  only  membrane  effects,  and  they  were  not  compared  with  any 
experimental  results.  The  analysis  by  Hu  et  al  was  an  inextensional 
analysis;  i.e.,  one  in  which  the  membrane  strains  are  identically  zero, 
with  solution  by  the  Rayleigh-Ritz  technique.  Their  results  appeared 
to  agree  with  their  previously  reported  experimental  results  (Reference 
18)  for  circumferential  wave  numbers,  n,  ranging  from  2  to  at  least  9. 
Sewall's  analysis  was  said  to  have  been  solved  by  a  Rayleigh-Ritz  tech¬ 
nique,  and  results  were  presented  for  the  lowest  unsynmetric  mode  only. 
These  results  agreed  quite  well  with  experimental  results  reported  by 
Mixson  (Reference  17).  Krause's  analysis  included  both  membrane  and 
bending  effects,  and  solution  was  accomplished  by  a  modified  Galerkin 
method  in  which  it  was  not  necessary  that  the  assumed  modal  functions 
satisfy  the  force  and  moment  boundary  conditions.  However,  the  results 
obtained  by  Krause  did  not  appear  to  give  better  agreement  with  the 
results  of  Reference  18  than  did  the  much  more  simplified  analysis  of 
Reference  42 . 

After  a  careful  review  of  the  analyses  mentioned  above,  it  was  decided 
to  use  the  Rayleigh-Ritz  method  of  analysis,  but  to  include  the  inexten- 
sional  modes  only.  This  choice  was  motivated  by  these  practical  consi¬ 
derations:  (1)  the  best  agreement  reported  to  date  for  the  unsymmetric 
vibrations  of  free-free,  homogeneous,  isotropic  shells  was  obtained  by 
the  inextensional  analysis  of  Hu  et  al  (Reference  42);  (2)  in  a  sandwich 
shell,  the  magnitude  of  the  membrane  strain  energy  relative  to  the 
bending  strain  energy  is  considerably  smaller  than  that  for  a  homogeneous 
shell.  (This  is  due  to  the  difference  in  the  ratio  of  membrane  stiffness 
to  flexural  stiffness.)  Thus,  in  a  sandwich  shell,  the  error  introduced 
by  neglecting  membrane  effects  is  small  except  for  the  axisymmetric  case 
(n  =  0) .  This  is  clearly  shown  by  an  approximate  relation  presented  by 
Platus  (Reference  44)  for  the  natural  frequencies,  w,  of  a  shell: 

a)  =  (w2  +  u?)1/2  (1) 

where  we  is  the  frequency  calculated  by  considering  only  extensional 
effects,  and  is  the  frequency  calculated  for  the  same  n  value  but 
considering  inextensional  effects  only.  In  deriving  Equation  (1) ,  it 
is  assumed  that  the  modal  functions  associated  with  the  two  kinds  of 
modes  are  geometrically  similar,  so  that  they  affect  the  frequency  only 
through  their  effect  on  the  stiffness.  The  hypothesis  that  the  exten¬ 
sional  effects  are  relatively  small  for  a  sandwich  conical  shell  were 
verified,  for  the  case  of  freely  supported  edges,  by  the  analysis  of 
Bacon  and  Bert  (Reference  40) (see  Figure  4). 

The  analysis  presented  in  the  next  section  considers  sandwich  shells 
in  the  foim  of  truncated  cones  with  orthotropic  facings  and  a  perfectly 
rigid  core.  All  components  of  translational  inertia  are  included,  but 
rotatory  inertia  is  neglected. 
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Figure  4.  Effects  of  Inextensional  and  Extensional  Deformations  on 
Combined  Frequencies  for  Free-Free  Conical  Shells  of 
Homogeneous  and  Sandwich  Types. 
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More  refined  analyses  are  currently  in  progress  at  the  University  of 
Oklahoma.  These  include  memhrane  and  bending  effects  in  the  ortho  tropic 
facings,  transverse  shear  effects  in  the  core,  and  all  components  of 
translational  and  rotatory  inertia.  Siu  (Reference  45)  is  considering 
free-free  boundary  conditions  and  is  using  the  energy  method  of  solu¬ 
tion.  Wilkins  (Reference  46)  is  using  the  Galerkin  method  and  is 
including  the  following  boundary  conditions  at  both  ends:  freely  sup¬ 
ported,  clamped,  and  free. 

b.  Analysis  of  Orthotropic  Sandwich  Conical  Shell  With  Free  Edges 

Hu,  Gormley,  and  Lindholm  (Reference  42)  presented  an  inextensional 
analysis  of  a  homogeneous,  isotropic,  conical  shell  with  free  edges. 
Following  the  inextensional  concept  pioneered  by  Lord  Rayleigh  (Reference 
42),  they  set  the  middle- surface  strain  components  equal  to  zero  and 
solved  the  resulting  set  of  first-order  differential  equations  to  obtain 
the  following  expressions  for  the  inextensional  displacements: 

u  =  A  sin  a  cos  a  sin  n0  cos  ai^t 

v  =  (A  +  Bs/s2)  n  cos  a  cos  n0  cos  u>^ t  (2) 

w  =  [A(n2  -  sin2  a)  +  Bn2  s/s2]  sin  n0  cos  u^t. 

Here  the  same  general  approach  as  used  by  Hu  et  al  is  employed,  except 
that  an  orthotropic  sandwich  conical  shell  is  considered  (see  Figure  5). 
Since  the  extensional  contribution  to  the  strain  energy  is  much  smaller 
relative  to  the  inextensional  controbution  in  a  sandwich  shell  than  in 
a  homogeneous  one,  an  inextensional  analysis  is  even  more  appropriate 
here  than  for  a  homogeneous  shell  (see  Figure  4). 

The  homogeneous,  linear  algebraic  equations  for  the  constants  A  and  B 
appearing  in  Equation  (2)  obtained  from  application  of  the  Rayleigh-Ritz 
technique  are  as  follows: 

(>2  rll  -  dU>  A  *  C»2  c12  -  d12>  B  -  0  (J) 

C12  "  d12^  A  +  C22  ’  d22^  B  =  0 
The  eigenvalues,  X,  for  these  equations  are  defined  as  follows: 

X  =  a,i  R^  (ph) 1/2 [ (1  +  n'2  cos2  a)/(n2  -  1)D0]1/2  (4) 

The  solutions  of  Equations  (3)  are 
X  =  [Q  ±  (Q2  -  4PR) 1/2]/2P 
and  the  corresponding  ratios  of  B  to  A  are 
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m  -  (a2  cu  -  djp/c*2  c12  -  a12) 


(5) 


where 


P  =  cn  c22  ■  c12 

Q  =  C11  d22  +  c22  dll  "  2c12  d12 
R  =  dll  d22  ‘  d12 

The  values  for  the  factors  in  Equations  (6)  are 

cn  =  [(1  -  y2)/2]  (1  -  [  (2n2  -  1)  sin2  a]/[n2(n2  +  cos2  a)]} 
c12  =  [(1  -  y3)/3][1  -  (sin2  c0/(n2  +  cos2  a)] 
c22  =  (1  -  y4)/4 

dn  =  {[(1/y2)  -  13/2H1  +  (DsQ/De)[(sin  a)/n2] 

=  (y'1  -  1) 


d22  =  loge  (1/y) 


where  y  is  the  completeness  parameter  (R-./R-) ,  and  for  a  sandwich  with 
thin  facings, 

De  -  Ee  tfh2/[2(l  -  vse  v0s)]  (8) 


Dse  '  “seV2 

Equations  (7)  are  identical  to  those  given  by  Hu  et  al,  except  the 
expression  for  d,, .  However,  for  the  isotropic  case  considered  by  Hu 
et  al,  D  n/D.  becomes  2(l-v),  and  the  expression  for  d, ,  coincides  with 
that  givift  by  Hu  et  al.  11 

When  the  expression  for  the  B/A  ratio,  Equation  (5) ,  is  substituted  into 
the  expression  for  the  circumferential  displacement,  the  second  of  Equa¬ 
tions  (2),  the  resulting  expression  for  the  normalized  modal  shape  is 

V  =  (1  ♦  [(X2  cn  -  dn)/(X2  c12  -  d12)](s/s2)}  (10) 

Similarly,  when  Equation  (5)  is  substituted  into  the  expression  for 
normal  displacement,  the  third  of  Equations  (2),  the  resulting  expression 
for  the  normalized  modal  shape  is 
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W  =  n2  U  -  (sin2  a)/n2  +  [(X2  cn  -  dn)/(A2  cu  -  d12)]s/s2) 

(ID 

Since  the  middle  surface  extensional  strains  and  the  meridional  bending 
strain  are  zero,  the  only  strain  existing  is  the  circumferential  bending 
strain,  e^Q.  The  circumferential  surface  bending  strain  is 

eb0  =  [fh  +  V/21  Ke  (12) 

where  is  the  change  in  circumferential  curvature. 

Thus,  the  meridional  distribution  of  the  dimensionless  circumferential 
surface  bending  strain  is 

eb0  "  fCSi/s)2  +  (B/A)(s2/s2s)]  (13) 

Although  the  results  of  this  simple  analysis  are  in  algebraic  form,  for 
the  values  of  geometric  and  material  parameters  involved,  subtraction 
of  two  large  and  nearly  equal  numbers  is  involved.  Therefore,  to  obtain 
reasonable  numerical  results,  the  above  equations  for  the  eigenvalues, 
eigenvectors,  modal  shape,  and  strain  distribution  were  programmed  in 
FORTRAN  IV,  double  precision  language  and  were  computed  on  the  Univer¬ 
sity  of  Oklahoma's  IBM  360/40  digital  computer.  The  program  was  written 
to  obtain  the  various  quantities  as  a  function  of  the  circumferential 
wave  number,  n,  with  the  semi-vertex  cone  angle,  a,  as  a  variable. 
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6.  EXPERIMENTAL  AND  ANALYTICAL  RESULTS 

The  experimentally  measured  values  of  resonant  frequencies  and  associated 
strain  distributions  and  damping  are  presented  in  this  section.  The  reso¬ 
nant  frequencies  and  associated  strain  distributions  for  the  two  lowest 
unsymmetric  modes  are  compared  with  the  analyses  presented  in  Section  5. 

Resonant  frequencies  for  various  meridional  modes  are  shown  in  Figure  6  as 
a  function  of  circumferential  mode  nunber.  On  this  figure,  the  two  curves 
that  were  derived  from  the  inextensional  analysis  presented  in  Section  5 
are  shown  along  with  the  experimental  resonant  frequency  points.  Each 
resonant  point  was  obtained  from  a  modified  Kennedy -Pancu  plot.  Some  typi¬ 
cal  plots  that  were  used  to  separate  the  resonant  frequencies  are  shown  in 
Figures  7  through  14.  As  can  be  seen  from  the  figures,  the  separation  of 
the  resonant  frequencies  that  were  very  close  together  was  accomplished 
quite  adequately  by  means  of  the  Kennedy- Pancu  technique.  Damping  factors 
were  determined  from  these  figures  using  the  technique  outlined  in  the 
appendix. 

It  can  be  seen  in  Figure  6  that  the  experimentally  measured  resonant- 
frequency  points  associated  with  the  two  lowest  unsymmetric  modes  agree 
quite  closely  with  the  frequencies  calculated  by  the  inextensional  analysis 
over  the  entire  range  of  circumferential  wave  numbers  covered. 

The  resonant- frequency  curves  for  the  higher  modes  indicate  the  same  trend 
as  those  of  the  two  lowest  modes.  In  the  region  of  low  circunferential  wave 
numbers,  n,  the  higher  modes  were  difficult  to  detect  because  the  resonant 
frequency  coincided  with  the  resonant  frequency  of  another  mode.  At  these 
lower  resonant  frequencies,  the  response  of  the  mode  having  a  higher  n 
would  dominate  that  of  the  lower-n  mode.  This  problem  was  one  of  the  major 
difficulties  encountered  in  applying  the  Kennedy- Pancu  technique  in  this 
investigation.  It  was  beyond  the  scope  of  this  investigation  because  of  the 
limitation  of  equipment,  but  a  solution  to  this  problem  would  be  to  intro*- 
duce  additional  excitation  systems.  By  placing  the  exciters  near  a  point 
that  is  a  node  of  the  strong  mode,  the  system  would  exhibit  a  larger  response 
in  the  weak  mode. 

Also,  the  higher  modes  displayed  the  same  trend  as  the  lower  modes  at  the 
higher  circumferential  wave  numbers  and  indicated  a  possible  convergence  of 
all  the  modes,  if  higher  wave  numbers  and  frequencies  could  have  been  ob¬ 
tained. 

At  the  lower  wave  numbers,  it  is  believed  that  the  higher  modes  will  not 
indicate  as  severe  of  a  hook,  similar  to  the  curves  presented  by  Watkins 
and  Clary  (Reference  47) ,  Hu  et  al  (Reference  18) ,  and  others  in  reporting 
the  unsymmetrical  resonant  frequencies  of  homogeneous -material  conical  shells. 
ITie  reason  for  this  was  discussed  in  Section  5  in  connection  with  Equation 
(1). 

Shown  as  Figures  15  and  16  are  the  plots  of  circumferential  strain  associated 
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Figure  15.  Normalized  Circumferential  Strain  Distribution  Along 
Meridian,  Inextensional  Modes,  n  =  2  through  5. 
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Figure  16.  Normalized  Circumferential  Strain  Distribution  Along 
Meridian,  Inextensional  Modes,  n  =  6  through  9. 
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with  the  two  lowest  meridional  modes.  Also  included  as  Figures  15  and  16 
are  the  analytical,  circumferential  strain  distributions  derived  from  the 
inextensional  analysis  presented  in  Section  5.  In  general,  there  is  very 
good  agreement  between  the  analytical  and  experimental  strain  distributions 
for  these  modes.  The  effect  of  exciter  mass  and  location  is  noticeable  in 
the  meridional  modal  shapes.  However,  this  is  to  be  expected  since  it  was 
observed  by  Mixson  (Reference  17) .  The  circumferential  modal  shapes  are 
symmetrical  about  the  exciter  location.  Moving  the  exciter  to  a  new  circum¬ 
ferential  location  on  the  shell  resulted  in  the  same  modal  shape  but  with 
the  nodal  points  shifted  circumferentially. 

Figure  17  shows  a  plot  of  the  meridional  strain  distribution  associated  with 
the  two  lowest  meridional  modes.  The  meridional  strain  distribution  is 
identically  zero  in  the  inextensional  analysis  presented  in  Section  5.  Even 
though  the  curves  indicate  a  meridional  strain,  the  values  are  small  and 
can  be  attributed  to  the  shell's  being  restrained  from  completely  free  move¬ 
ment  at  the  shaker- attachment  location  and  at  the  suspension  points. 

Since  the  experimental  data  and  the  values  derived  by  the  inextensional 
analysis  are  in  close  agreement,  it  is  assumed  that  the  lowest  experimental 
unsymmetrical  modes  are  the  inextensional  modes.  The  agreement  is  close 
for: 

(a)  the  resonant  frequencies  as  a  function  of  circumferential  wave 
number, 

(b)  the  circumferential  strain  as  a  function  of  meridional  position, 

(c)  the  meridional  strain  as  a  function  of  meridional  position. 

Hereafter,  these  two  lowest  experimental  mode'  will  be  referred  to  as  the 
inextensional  modes,  since  these  agree  closrly  with  the  analytical  inexten¬ 
sional  modes.  The  other  modes  which  have  higher  resonant  frequencies  will 
be  referred  to  as  general  modes,  since  they  exhibit  both  extensional  and 
inextensional  deformation  as  discussed  in  Section  5. 

The  circumferential  strain  distribution  patterns  associated  with  the  circum¬ 
ferential  modes  for  the  inextensional  and  general  modes  are  shown  in  Figures 
18  and  19.  The  circumferential  strain  distributions  associated  with  the 
general  modes  are  shown  in  Figure  20.  The  distributions  of  meridional 
strain  along  the  meridian  are  shown  in  Figure  21.  It  can  be  noted  from 
this  figure  that  the  meridional  strain  is  not  uniform  along  the  meridian, 
in  contrast  to  the  cases  of  the  two  inextensional  modes.  It  can  be  seen 
in  Figure  6  that  the  resonant  frequency  for  the  m  =  1  general  mode  is  close 
to  inextensional  mode  B  at  the  low  circumferential  wave  numbers,  but  as  the 
wave  number  increases,  the  curves  separate  and  the  m  =  1  mode  converges  to 
the  m  =  2  general  mode.  Also,  it  can  be  noted  in  this  figure  that  all  of 
the  general  modes  tend  to  converge  at  the  higher  frequencies.  This  conver¬ 
gence  is  characteristic  of  the  dynamics  of  both  cylindrical  and  conical 
shells. 
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Figure  17.  Meridional  Strain  Distribution  Along  Meridian,  Inextensional 
Modes,  n  =  2.  3.  6.  8. 
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Figure  20.  Normalized  Circumferential  Strain  Distribution  Along  Meridian, 
General  Modes. 
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The  resonant- frequency  data  and  the  associated  strain  distributions  indi¬ 
cated  that  sane  coupling  between  modes  occurred.  These  coupled  modes  were 
a  combination  of  two  circumferential  modes  having  the  same  meridional  mode 
when  the  resonant  frequencies  were  close. 

Pronounced  coupling  phenomena  were  observed  at  resonant  frequencies  of  33 
cps,  64  cps,  and  94  cps.  In  all  of  these  cases,  the  strain  distributions 
were  typically  the  same  except  for  the  difference  in  circumferential  wave 
numbers.  Shown  as  Figures  22  through  24  are  typical  plots  of  the  nodal 
patterns  associated  with  the  above  frequencies.  In  each  of  these  figures, 
the  nodal  patterns  were  symmetrical  about  the  exciter  but  not  evenly  spaced. 

The  phenomenon  of  U-shaped  nodal  lines,  i.e.,  the  existence  of  more  circum¬ 
ferential  waves  at  the  large  end  of  a  free- free  conical- frustum  shell  than 
at  the  small  end,  was  observed  in  experiments  reported  by  Watkins  and  Clary 
(Reference  16).  This  has  resulted  in  considerable  controversy,  since  they 
were  not  observed  at  all  in  any  of  the  similar  experiments  reported  by  Hu 
et  al  (Reference  18) .  Hu  (Reference  48)  suggested  that  the  U-shaped  modes 
were  caused  by  shaker  effects,  since  his  experiments  (Reference  18)  were 
conducted  by  magnetic  fields  (no  direct  contact  with  shell) .  Additional 
comments  were  made  by  Koval  (References  49  and  50) ,  who  showed  photographs 
of  similar  phenomena  in  vibrational  experiments  on  cylindrical  shell.  He 
attributed  the  U-shaped- mode  phenomenon  to  the  dynamic  asymmetries  caused 
by  the  spot -welded  longitudinal  seams  present  in  Watkins  and  Clary’s  shells. 
According  to  the  theory  originated  by  Tobias  (Reference  51)  and  applied  by 
Arnold  and  Warburton  (Reference  52),  the  presence  of  small  imperfections 
causes  each  of  the  two  close  resonant  frequencies  to  have  a  preferential 
nodal  pattern.  Then,  under  excitation  at  a  frequency  between  the  two  com¬ 
ponent  frequencies,  the  two  different  component  modes  combine  to  produce 
the  complicated  mixed-mode  pattern. 

As  a  result  of  the  controversy,  following  the  suggestion  of  Watkins  and 
Clary  (Reference  47),  Mixson  (References  17  and  53)  conducted  an  extensive 
series  of  additional  experiments  on  shells  with  and  without  longitudinal 
seams.  In  these  experiments,  he  used  air  shakers,  as  well  as  electro- 
dynamic  ones,  and  varied  the  ratio  of  the  frequency  of  the  axial  rigid-body 
mode  (controlled  by  the  suspension  stiffness)  to  that  of  the  lowest  elastic 
shell  mode.  He  concluded  that  imperfections  were  probably  the  principal 
cause  of  the  mixed,  modes,  since  they  occurred  most  often  in  the  cones  with 
seams  and  smaller  wall  thickness  (in  which  it  is  more  difficult  to  control 
imperfections) . 

The  coupling  observed  was  probably  caused  by  the  presence  of  the  four 
equally  spaced  longitudinal  lap  joints  on  the  inside  and  outside  facing 
(i.e.,  a  total  of  eight  lap  joints;  one  every  45  degrees).  However,  it  is 
interesting  to  note  that  the  only  modes  which  coupled  were  inextensional 
modes  (A  and  B).  There  was  no  observed  coupling  of  the  general  modes 
(m  *  1,  2,  3)  with  each  other  or  with  the  inextensional  modes.  This  may 
have  been  because,  for  the  particular  shell  construction  and  geometry,  the 
inextensional  modes  were  the  only  ones  sufficient^  close  together  (in 
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frequency  as  well  as  in  circumferential  wave  nunber)  to  couple. 

At  present,  the  authors  know  of  no  analytically  based  criterion  for  deter¬ 
mining  a  priori  whether  or  not  two  adjacent  modes  will  couple  to  form  a 
combined,  mode.  Since  the  available  experimental  data  are  quite  limited  and 
expensive  to  obtain,  an  analytically  based  criterion,  rather  than  an  empiri¬ 
cally  based  one,  would  be  quite  useful. 

Since  the  conical  shell  used  here  had  a  small  cone  angle,  it  could  be 
approximated,  with  reasonable  accuracy,  by  a  cylinder  for  calculating  the 
frequencies  of  certain  higher  modes  peculiar  to  sandwich  structures.  For 
these  calculations  Yu's  analysis  (Reference  22)  was  used.  It  was  found 
that  the  calculated  frequencies  for  all  three  of  these  modes  (shear  modes 
in  both  the  meridional  and  circumferential  planes  and  the  thickness -normal 
mode)  were  much  higher  than  the  highest  frequencies  attained  in  the  experi¬ 
ments.  Therefore,  it  is  presumed  that  no  modes  of  these  types  were  excited 
in  the  experiments  reported. 

Figure  25  shows  a  plot  of  the  damping  logarithmic  decrement  as  a  function 
of  resonant  frequency.  The  values  for  this  figure  were  obtained  from  the 
Kennedy-Pancu  plots  and  converted  to  logarithmic  decrements.  As  pointed 
out  in  the  appendix,  the  damping  coefficient  of  the  system  can  be  obtained 
from  the  Kennedy-Pancu  plots.  The  damping  coefficient  is  converted  to 
logarithmic  decrement  by  multiplying  the  coefficient  by  2tt.  The  tendency 
for  the  damping  to  decrease  as  the  frequency  increases  agrees  with  data 
presented  by  Shoua  (Reference  54) .  The  damping  logarithmic  decrement  at 
the  higher  frequencies  agrees  with  values  presented  by  Shoua  and  Bert  et 
al  Reference  8)  for  beams  of  the  same  type  of  construction  and  materials 
as  used  in  the  present  shells  (sandwich  with  fiber  glass  epoxy  facings  and 
aluninum  honeycomb  core) . 

The  same  damping  logarithmic  decrement  is  shewn  in  Figure  26  as  a  function 
of  circumferential  wave  number. 
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Figure  25  .  Damping  Logarithmic  Decrement  as  a  Function  of  Resonant 
Frequency. 
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Figure  26.  Damping  Logarithmic  Decrement  as  a  Function  of  Circumferential 
Wave  Number. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS 


a.  Conclusions 

(1)  With  sufficiently  accurate  electronic  instrumentation  now  avail¬ 
able  copmercially  and  with  sufficient  care  in  conducting  the 
experiments ,  a  modified  Kennedy-Pancu  technique  has  been  applied 
successfully  to  shell- type  structures  to  determine  very  accurately 
the  resonant  frequencies,  modal  strain  distributions,  and  damping 
associated  with  pure  modes. 

(2)  For  free-free  sandwich  shells,  an  inextens ional  analysis  predicts 
quite  accurately  the  resonant  frequencies  associated  with  the 
first  two  unsymmetric  modes.  Under  the  conditions  of  the  present 
experiments,  extens ional, core  shear,  and  core  normal  effects  can 
be  neglected  for  these  two  modes. 

(3)  The  combined  modes  with  U-shaped  nodal  patterns,  as  reported  in 
previous  experiments  on  free- free  cones,  were  observed.  However, 
in  ail  cases,  they  represented  coupling  between  inextensional  modes 
only. 

b.  Recommendations 

(1)  For  applications  involving  very  high  frequency  excitations,  such 
as  blast  loading  or  high-frequency  sonic  excitation,  the  higher 
order  effects  of  transverse  shear  deformation,  thickness  normal 
deformation,  and  rotatory  inertia  may  be  important.  To  investigate 
higher  frequencies,  additional  research,  using  the  same  size  of 
sandwich  shells  and  larger  excitation  force,  would  be  required. 

(2)  An  analysis  is  needed  to  predict  a  priori  whether  or  not  two  adja¬ 
cent  modes  will  or  will  not  couple  into  combined  modes. 
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APPENDIX 


MODIFIED  KENNEDY -PANCU  TECHNIQUE 

Kennedy  and  Pancu  (Reference  55)  introduced  their  method  to  obtain  the 
uncoupled  resonant  frequencies  and  associated  modal  shapes  in  1947  in  con¬ 
nection  with  ground  vibration  tests  on  aircraft.  Their  proposals  were  made 
entirely  on  the  basis  of  theoretical  arguments,  with  experimental  results 
mentioned  but  not  reported  in  detail.  Since  the  appearance  of  their  paper, 
the  technique  has  been  used,  but  only  a  limited  number  of  papers  on  the 
technique  have  appeared.  All  the  presentations  on  the  technique  were 
oriented  toward  beam- type  structures,  none  to  the  authors'  knowledge  deal¬ 
ing  with  shell- type  structures.  The  technique  was  used  with  success  in 
airplane  ground  vibration  tests  (Reference  55),  with  airplane  inflight 
flutter  tests  (Reference  56),  and  with  beams  (Reference  57).  The  present 
application  is  the  first  known  use  for  shell  structures.  The  technique  is 
especially  pertinent  to  shell  structures,  etc.,  because  of  the  closely 
coupled  resonant  frequencies  that  are  often  encountered  in  shells. 

The  main  limitation  in  using  the  technique  in  vibration  testing  is  the 
accuracy  with  which  the  phase  between  the  excitation  force  and  the  output 
force  is  measured.  When  Kennedy  and  Pancu  published  their  paper,  there  was 
no  electronic  equipment  in  existence  that  could  measure  this  phase  relation 
to  the  necessary  accuracy.  Even  with  accurate  phase-measuring  devices, 
there  are  a  number  of  limitations  in  using  the  Kennedy- Pancu  technique. 

As  pointed  out  by  Pendered  and  Bishop  (Reference  58)  and  also  by  Bishop 
and  Gladwell  (Reference  56),  near  resonance  there  is  a  rapid  change  in 
phase  with  respect  to  change  in  frequency,  making  the  measurement  of  this 
parameter  difficult.  Since  this  is  a  very  important  parameter  in  deter¬ 
mining  the  resonant  point,  its  measurement  must  be  accurate.  In  the  present 
work,  this  difficulty  was  overcome  by  developing  a  testing  procedure  in 
which  the  phase  variation  was  used  as  the  control  parameter.  A  highly 
accurate  digital  phase  meter  was  used  to  measure  the  phase. 

Distortion  in  the  output  signals  causes  errors  in  the  off-resonant  frequen¬ 
cies.  When  the  signal- to-noise  ratio  becomes  small,  the  phase  signal  becomes 
distorted  with  undesired  signals.  Fortunately,  this  does  not  present  too 
much  difficulty,  since  the  main  area  of  interest  is  near  resonance  where  the 
signal- to-noise  ratio  is  the  largest.  This  should  be  enough  to  determine 
the  characteristic  circular  arc  to  deduce  the  resonant  frequency  and  the 
damping.  The  Lissajous  pattern  was  used  to  check  the  quadrant  of  the  phase 
enabling  the  determination  of  the  circular  arc  extremities. 

The  Kennedy-Pancu  analysis  was  derived  for  a  system  that  was  linear.  In 
nearly  all  practical  systems,  some  nonlinearities  are  inherent.  This  was 
minimized  in  the  present  tests  by  establishing  an  output  level  and  varying 
the  input  force  amplitude.  In  using  the  Kennedy-Pancu  method,  the  Argand 
plots  (amplitude  versus  phase)  are  normalized  by  the  input  force  level;  i.e., 
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complex  receptance  is  plotted  on  the  complex  plane.  Thus,  resonance  occurs 
when  the  output  force  is  the  largest  or  the  input  force  is  the  smallest. 
Holding  the  output  level  constant  established  the  same  displacement  on  the 
shell  throughout  the  frequency  range,  thus  minimizing  the  nonlinearities, 
regardless  of  whether  they  are  geometric  (large  deflections)  or  material 
(nonlinear  stress-strain  and  nonlinear  damping) . 

The  Kennedy-Pancu  technique  has  an  advantage  over  other  methods  for  deter¬ 
mining  resonant  frequencies  because  it  is  able  to  separate  modes  that  are 
closer  to  each  other  and  also  able  to  come  closer  to  determining  the  pure 
modal  shape  (i.e.,  without  coupling  with  other  modes). 

In  applying  their  technique  to  a  multi-degree-of-freedom  system  (as  appli¬ 
cable  to  this  case) ,  Kennedy  and  Pancu  made  the  assumption  that  each  mode 
acts  independently  of  the  others;  i.e.,  each  mode  responds  to  the  applied 
force  like  a  single- degree- of- freedom  system.  Also,  they  assumed  that  the 
system  is  linear  and  that  superposition  of  modes  will  hold.  In  other  words, 
they  assumed  the  existence  of  normal  modes,  which  is  a  widely  accepted 
assumption  in  multi-degree-of-freedom  vibrational  analysis. 

The  response  of  a  damped  single- degree -of- freedom  system,  with  Kimball- 
Lovell-type  material  damping  and  excited  by  a  simple  harmonic  force  of 
amplitude  F  and  frequency  w,  can  be  expressed  as  follows 


X  =  [(k  -  mw2)2  +  b2] 

(14) 

tan  <J>  =  b/ (k  -  mu2) 

(15) 

where 

b  =  material  damping  factor* 
k  =  spring  rate 
m  =  mass 
X  =  X0/F 

x  =  displacement,  given  by  Xq  sin  (wt  -  $) 

<{>  =  phase  angle  between  the  displacement  and  forcing  function 

If  a  complex-plane  plot  is  made  of  the  complex  receptance,  Xe1'*’,  with  the 
mass,  m,  damping  factor,  b,  and  spring  rate,  k,  fixed  and  the  frequency 
varied,  the  result  is  as  shown  in  Figure  27. 


*The  material  damping  factor  is  not  to  be  confused  with  the  viscous  damping 
coefficient,  which  is  equal  to  b/w. 
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Figure  27.  Argand  Plot  for  Single-Degree-of-Freedom  Damped  System. 


The  point  where  the  curve  crosses  the  imaginary  axis  is  the  point  in  which 

w2  =  w2  =  k/m 

and  is  defined  as  the  resonant  frequency  of  the  system.  Then  at  the  point 
of  resonance,  the  distance  from  the  origin  to  the  intercept  of  the  imaginary 
axis  is  related  to  the  material  damping  factor,  b,  as  follows: 

X  =  F/b 
o 

Kennedy  and  Pancu  noted  that  the  Argand  plot  traced  a  circular  arc  near  the 
resonant  point  and  that  a  circle  could  be  passed  through  the  points.  By 
drawing  the  circle,  the  resonant  frequency  of  the  system  would  occur  at  the 
point  in  which  the  phase  angle  made  90  degrees  with  the  forcing  function  and 
also  the  line  would  pass  through  the  center  of  the  circle.  It  was  also 
noted  that  the  resonant  point  occurred  when  the  rate  of  change  in  frequency 
with  respect  to  phase  angle  was  the  greatest. 

The  procedure  for  obtaining  the  damping  coefficient  from  the  Argand  plot 
can  be  simplified  by  eliminating  the  necessity  of  making  absolute  calibra¬ 
tions  of  the  excitation  and  response  transducers.  This  is  accomplished  by 
expressing  the  damping  coefficient,  b,  in  terms  of  frequency  as  follows 

b/k  =  (oij  -  u)2)/u)n  (16) 

where  b/k  is  a  dimensionless  damping  factor  and  u-  and  are  the  frequencies 
at  the  half -power  points  on  the  response  curve.  Further  assuming  that  the 
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response  curve  is  symmetrical  about  the  resonant  frequency,  u>  ,  then 
Equation  (16)  becomes  n 

b  =  2(1  -  (17) 

The  half -power-point  frequency  can  readily  be  determined  from  the  Argand 
plots  by  drawing  a  line  from  the  center  of  the  circle  along  a  45-degree 
line;  the  point  of  intersection  of  this  line  with  the  circle  is  the  half¬ 
power-point  frequency,  cu..  Applying  these  values  to  Equation  (17),  the 
results  will  give  the  uncoupled  damping  coefficient  for  the  system. 

When  more  than  one  resonance  is  present  in  the  vibrating  system,  then  each 
will  introduce  an  influence  on  the  response  data.  The  closer  the  resonant 
frequencies  are  together,  the  greater  will  be  the  influence. 

Using  the  assumption  that  each  mode  acts  independently,  then  the  system 
can  be  defined  by  a  series  of  Equations  (14)  and  (15),  and  the  Argand  plot 
of  the  system  will  be  a  series  of  circles,  each  circle  crossing  its  own 
imaginery  axis  at  its  resonant  frequency.  The  plot  will  be  askew  depending 
on  how  close  together  the  resonances  are,  but  each  resonance  will  act  inde¬ 
pendently  and  will  have  its  own  circle  on  a  displaced  set  of  coordinate 
axes.  By  applying  the  same  reduction  to  each  circle,  the  system  can  be 
separated  into  individual  single-degree-of-freedom  systems  acting  indepen¬ 
dently.  It  should  be  pointed  out  that  for  a  continuous  system  the  Kennedy - 
Pancu  technique  is  applicable  because,  even  though  the  system  is  continuous, 
the  resonant  frequencies  are  discrete. 

Once  the  Argand  plot  is  made,  circular  arcs  are  drawn  through  the  data  points 
that  indicate  a  relatively  large  change  in  frequency  with  relation  to  change 
in  phase.  By  drawing  a  line  that  is  displaced  90  degrees  from  the  forcing 
function  and  passes  through  the  center  of  the  circle,  the  intercept  of  this 
line  with  the  circle  determines  the  resonant  frequency  and  the  diameter  (XQ) 
of  the  circle  are  functions  of  the  damping  factor,  as  mentioned  before. 

Each  circle  on  the  plot  indicates  the  resonant  frequency  and  damping  of  the 
individual  mode  acting  alone  can  be  obtained  from  the  distance  from  the 
center  to  the  edge  of  the  circle  for  the  individual  frequency  points. 
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